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Abstract 

We question the validity of the grand canonical ensemble for the description 
of Bose-Einstein condensation of smah ideal Bose gas samples in isolated har- 
monic traps. While the ground state fraction and the specific heat capacity 
can be well approximated with the help of the conventional grand canonical 
arguments, the calculation of the fluctuation of the number of particles con- 
tained in the condensate requires a microcanonical approach. Resorting to the 
theory of restricted partitions of integer numbers, we present analytical and 
numerical results for such fluctuations in one- and three-dimensional traps, 
and show that their magnitude is essentially independent of the total particle 
number. 
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I. BOSE EINSTEIN CONDENSATION IN ISOLATED TRAPS: CAN WE RELY 
ON THE GRAND CANONICAL ENSEMBLE? 



After the breakthrough in the preparation and detection of Bose-Einstein condensates 
of rubidium sodium [Q, and hthium the second generation of experiments on ultra- 
cold, magnetically trapped Bose gases now starts to probe condensate properties such as 
the temperature dependence of the condensate fraction or collective excitations of the 
condensate A matter-of-principle experiment has shown the feasibility of an output 

coupler for condensed trapped atoms |P, two different, overlapping condensates have been 

1^], and interference fringes between two freely expanding con- 
thus proving first-order coherence of "mesoscopic" 
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produced in a single trap 
densates have been demonstrated 
matter waves. Higher-order coherence has been inferred from three-body recombination 
rates of condensed atoms |0, and even the propagation of sound in a Bose condensate 



could be directly observed |jT3| . 



These developments naturally force the theorist to investigate a conceptually important 
question: When dealing with Bose-Einstein condensates consisting of about 10^ to 10^ 
particles in a magnetic trap, to which extent can the different thermodynamical ensembles 
be considered as equivalent? A Bose gas cooled below the condensation point and kept in 
a trap is neither in thermal contact with a heat bath, nor exchanging particles with some 
particle reservoir. Hence, the only ensemble that is appropriate for describing this situation 
is the microcanonical one. And yet, the theoretical discussion of trapped condensates so far 
mainly relies on the grand canonical ensemble. 

Within the grand canonical ensemble, the description of a trapped, ideal Bose gas con- 
sisting of particles is comparatively simple. We consider an isotropic harmonic trap with 
oscillator frequency u, so that the degree of degeneracy of the j-th single-particle state is 



gj = -{j + i){j + 2). 



(1) 



Hence, the difference /i of the chemical potential and the single-particle ground state energy 
is determined by the equation 



J 72 + 3j/2 + 1 



(2) 



where ks denotes the Boltzmann constant, and T is the temperature. This sum can be 
evaluated approximately by introducing the density of states 



PiE) 
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with 7 = 3/2, which follows immediately from eq. ([1|). Converting the sum into an integral 



and applying standard arguments, one then finds the condensation temperature |14-1(; 
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where 
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(5) 



denotes the condensation temperature pertaining to the large- A^-hmit [|T^,|I^; (^{z) is the 
Riemann zeta function. The lowering of T^^ with respect to Tq^^ is of order N~^^^ and 
results from the enhancement of the density of states (|^) above the leading "volume" -term: 
since there are more states available, the need to condense arises only at lower temperatures. 
It should be noted that, strictly speaking, there is no well-defined condensation temperature 
for a gas consisting of a finite number of particles. However, Fig. 1 clearly demonstrates 
that the onset of condensation in a three-dimensional harmonic trap becomes quite sharp 
already for particle numbers of the order of 10^, and eq. (1) describes this onset very well. 

In contrast to the textbook case of the free Bose gas, the heat capacity of the harmonically 
trapped Bose gas exhibits a steep drop at the condensation point. As shown in Fig. 2, this 
drop becomes a discontinuous jump by about 6.6 ks per particle in the large- A^-limit. 

Remaining within the scope of the grand canonical ensemble, this analysis can be made 



formally more precise, and generalized to moderately anisotropic harmonic traps |T^-p2|, 
which changes the value of 7. But can one rely on the grand canonical ensemble? Ziff, 
Uhlenbeck, and Kac have made a case that the grand canonical ensemble does not represent 
any physical situation for the condensed ideal Bose gas, and advocate that its anomalous 
predictions should be ignored These "anomalous predictions" are related to what we 



call the "grand canonical fluctuation catastrophe": within the grand canonical ensemble, 
the r.m.s. -fluctuations 6Nq of the ground state occupation number Nq are given by [p4| , ^ 



i6Nof = No (No + 1) , (6) 

implying that 6N0 is of the order of the total particle number N below the condensation 
point. However, if we consider an isolated Bose gas in a trap, all particles occupy the ground 
state at zero temperature, so that the true fluctuations of A^o vanish. The grand canonical 
prediction for 6N0 thus differs drastically from the microcanonical one. 



II. MICROCANONICAL APPROACH TO FLUCTUATIONS OF THE GROUND 

STATE OCCUPATION 

Since the magnitude of the fluctuations 6N0 is related to the coherence properties of the 
condensate, it is of substantial interest to compute the true, i.e., microcanonical fluctuations 
for a trapped ideal Bose gas. We will study an isotropic harmonic trapping potential with 
oscillator frequency u in d dimensions, and denote by n the number of excitation quanta for 
some preassigned value of the excitation energy E: 

The prime task now is to determine the number ^l^'^\n\N) of microstates. Since there are 
generally many microstates where only a part of the N particles carries all n excitation 
quanta, leaving the other particles in the ground state, ^l^'^\n\N) equals the number of 
possibilities for distributing the n quanta over at most N Bose particles. Then the difference 
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^l^^\n\Nex) — ^^'^^tT'INqx — 1) is the number of possibilities for distributing n quanta over 
exactly Nq^ particles, with A^^ex ranging from 1 to A^, so that 

is the probability for finding A'ex out of particles excited when the total excitation energy 
is n ■ hu. Since the remaining — A"ex particles occupy the ground state, the first moment 
(A'ex) of the distribution (|[) yields die microcanonical expectation value of the ground state 
occupation number according to 

(No) = N- (ATex) ; (9) 
the corresponding fluctuations follow from 

(^A^o)^ = {Nl) - (A^ex)^ . (10) 



III. CONDENSATE FLUCTUATIONS IN A ONE-DIMENSIONAL OSCILLATOR 

POTENTIAL 



The case d = 1 had been considered already in 1949 by Temperley p6[, and shortly 
thereafter by Nanda It is of particular interest, since it can be treated in detail analyti- 
cally with the tools furnished by the theory of the partion of integers: the number ^l^^\n\N) 
of microstates corresponds to the number of partitions of the integer n into at most A^ sum- 
mands; the commutativity of the summands reflects the indistinguishability of the Bosons. 
We now have to distinguish two cases: if n < A^, then the fact that the number of particles 
is finite has no consequences for the enumeration of microstates, and we are dealing with 
unrestricted partitions of n. If n > A^, the partitions of n are restricted by the requirement 
that the number of summands, corresponding to the number of excited particles, does not 
exceed A^. 

Let us first consider the case n < N. We follow the usual convention and denote the 
number of unrestricted partitions of n as p{n). Introducing the dimensionless inverse tem- 
perature ^ = huj/{kBT), it is an elementary exercise to show that 

oo oo -| 

gp(n)e-"^-n i_ (11) 



n=0 



which means that the generating function for p{n) corresponds, physically speaking, to 
the canonical partition function of a fictituous system of infinitely many distinguishable 
harmonic oscillators with frequencies that are integer multiples of u. Given this, it is natural 
to introduce the new variable x = exp(— ^) and to extract the numbers p{n) by inverting 
eq. (0) with the help of the saddle point approximation, after expressing ^ in terms of n. 
But we have to be careful: for temperatures that may physically be considered as "low", 
but are still high compared to hu/kB, the variable exp(— ^) is close to unity, so that we need 
to know the behaviour of Z^\$^) in the vicinity of a singularity. However, there exists the 
remarkable identity (see eq. (1.42) in ref. [E3]) 
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Z'^' (0 = / J exp - 1) Zg^ (4vr (12) 

that links the low-temperature behaviour of our fictituous oscillator system to its high- 
temperature dynamics. This allows us to derive the low-temperature approximation 

ln^^HO^^ + ^lne-^ln27r, (13) 

which, in turn, yields the desired energy-temperature relation: 

n = -§-^lr.Zg\0 - (14) 

Now we can compute the canonical entropy 

S{0/kB = n^ + \nZ2,\0 (15) 
and apply "Bethe's theorem" P3 to get 



, , exp[S{Un))/kB] 1 / FT] 

The right hand side is nothing but the celebrated Hardy-Ramanujan formula for the number 
of unrestricted partitions [^, so we have just sketched a physicist's solution to a number- 
theoretical problem. There is one detail that deserves particular attention: Bethe's theorem 



amounts to the inversion of the canonical partition sum (11) within the saddle point ap- 
proximation, so that the dominant contributions to this sum are properly taken into account 
within the Gaussian approximation. As a result, the microcanonical entropy Inp(n) differs 
from the canonical entropy S{^{n)) /kB by the saddle point correction — ln(— 27r9n/9^)/2. 
Such differences between thermodynamical quantities pertaining to different ensembles are 
characteristic for small systems. Expressed the other way round, the usually assumed equal- 
ity of thermodynamical quantites in different ensembles holds to the extent that such saddle 
point corrections can be neglected. 

Now we can turn to the case n > N, where the number of quanta exceeds the number 
of particles, and have to determine the number ^l^^\n\N) = pN{n) of restricted partitions 



of n. In principle, one can proceed as in the previous case, since there is the identity |^ 



N 



n=0 



exp(-jO 



The generating function for PN{n), that is, the canonical partition function for N ideal 
Bosons trapped by a one- dimensional harmonic potential, equals the canonical partition 
function of a system of N harmonic oscillators with frequencies u, 2uj, . . . , Nuj. The 
corresponding asymptotic formula for pN{n), which turns out to be rather intricate, has 
been given by Auluck and Kothari [Q. However, a beautiful theorem due to Erdos and 
Lehner helps to simplify the analysis: if, for some given n and x, the number obeys 
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^^v^nriT!^^^ with C'-Try'l, (18) 



Um = expf4e-4-) . (19) 



then 



p{n) \ C 

Hence, for given large values of n and we define x according to 



and obtain the approximation 



N Inn 

x=^- — , (20) 



P7v(^) exp(^-^e aC^j . (21) 

In order to check this approximation, we first determine the numbers PAr(n) = QS^\n\N) by 
means of the saddle point inversion of eq. (0), evaluate the distributions (|), and compute 
the relative fluctuations 5Nq/N according to eq. (0). The results for = 10^, 10^, and 
10^ are shown as full lines in Fig. 3. The corresponding data obtained with the help of 
the Erdos-Lehner approximation ( pi]) are indicated by the dashed lines. Evidently, the 
approximation is quite good already for = 10^. 

The merit of the approximation (0) lies in the fact that it allows us to derive an analytical 
expression for the microcanonical low-temperature ground state fluctuations of the ideal Bose 
gas trapped by a one- dimensional harmonic potential ||32||: 



As already seen in Fig. 3, the low-temperature fluctuations vanish linearly with temperature. 
Most notably, they do not depend on the particle number A^. This is immediately obvious 
for temperatures below the "restriction temperature" T^^"* that is defined by the condition 
n{Tj^^) = N, since below T^^^ the number of microstates becomes independent of A^, so that 
"the condensate has no chance to know how many particles it consists of" . However, for d = 1 
this restriction temperature is merely of the order of A^^/^, since k^Tj^^ ^ huj{N/({2)y^'^, 
wheras the A^-independence of SNq actually persists almost up to ksTo^^ = fiuj{N/ In N), 



where the occupation of the ground state becomes significant |^2|. This A^-independence 
of 6Nq appears to be characteristic for isolated condensates in general, as we shall indicate 
below. 



IV. CONDENSATE FLUCTUATIONS IN A THREE-DIMENSIONAL TRAP 



The analysis for the case d = 3 proceeds in close analogy to the case d = 1. The 
restriction temperature, determined from n{T^ ) = N, now reads 

Tj^^ ^0.8 N~^/^^tP , (23) 
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(3) 

which means that is about one fourth of the condensation temperature for a gas consist- 
ing of 10^ particles. Below T^^ the number of microstates, and hence all thermodynamical 
properties of the trapped gas, are independent of the particle number N. Accordingly, we 
write Q^^\n\N) = Q''^\n) for n corresponding to temperatures T less than so that 
fl^^\n) denotes the three-dimensional analogue of p{n). The generating function for fl^^\n) 
again equals the canonical partition function of a system of infinitely many, distinguishable 
harmonic oscillators PB[ |: 



n) e 



n=0 



n - 

fJl [1 - exp(-jO]* 



(24) 



The saddle point inversion of this generating function is straightforward, but a bit tedious. 



The result has been given by Nanda in 1951 [p^ : 



n 



-25/32 



,,,exp 4C(4)n3/4 + -C(3)ni/2 + 
47r(3C(4))'/' V 2^ 



^ 8 C(4) 



n 



(25) 



with 



n 



n 



3C(4) 



and 



B 



OSf _ C(2)C(3) 3 1 
8C(4)2 4C(4) 2^ ^ ^ ^ 2^ ^ ' 



This is one of the rare examples of a truly asymptotic formula for the number of microstates 
of a non-trivial Bose system. 

For temperatures above the finiteness of the particle number restricts the number of 
microstates, and we have to compute QS^\n\N). It turns out that the logarithm of Q!~^\n) 
provides a quite good approximation to ln(^fi'''^^(n|A^)^ even up to the condensation temper- 
ature |3^, so that the entropy of the fictituous Boltzmannian oscillator system described by 



the partition function (24) approximates the entropy of the trapped Bose gas for all tem- 
peratures below the onset of condensation. This means that the distributions pf^{N^^\n) 
introduced in eq. (||) must be well peaked; those microstates where the n > N quanta are 
actually "spread out" over all particles carry only a minor statistical weight below . 
But still, this finding is of no help for computing the fluctuations 6No, since this compu- 
tation requires, according to eq. (§), to consider auxiliary systems consisting of A^^ex < A^ 
particles and to determine the numbers Q^^'' {n\Ncx) ■ The problem now is that an immediate 
analogue to the generating function ( [T7| ) does not exist, that is, the canonical A^-particle 
partition function Z^\^) is not known in closed form. However, there exists the recursion 
formula |]35|-p^ 



N 



k=l 



id) 



(26) 
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that can numerically be evaluated and inverted within the saddle point approximation, so 
that we can determine the distributions p^{Nex\n) at least numerically with high preci- 
p3| . The results for = 1000 are shown in Fig. 4: as anticipated, the distributions 



sion 



are well peaked, and moreover remarkably close to Gaussians. 

From the expectation values of these distributions we obtain the microcanonical expec- 
tation values of the ground state occupation number, depicted as the full line in Fig. 5. 
The agreement with the corresponding grand canonical data (dashed) is stunning, but in 
view of the preceding Fig. 4 not unexpected: the predictions of the different ensembles will 
be the same as long as the most probable values of iVex coincide with the expectation val- 
ues, which is the CclSG clS long as the distributions p'^{Nex\n') stay symmetrical. This, in 
turn, is guaranteed as long as the support of these distributions stays away from the border 
Nex/N = 1. Hence, there is a slight difference between the microcanonical and the grand 
canonical ground state fraction only in the immediate vicinity of the condensation point. 
For d = 1 this difference is much more pronounced, since p^{NQx\n) is considerably less well 
peaked than its counterpart for (i = 3. 

Fig. 6 shows the microcanonical fluctuations of the ground state occupation number for 
d = 3 and N = 200, 500, and 1000. Remarkably, the low-temperature fluctuations for 
these three systems agree perfectly. More generally, if we compare two trapped Bose gases 
with different N under otherwise identical conditions, then their fluctuations 6Nq are nearly 
identical for temperatures below the lower of their condensation temperatures. This finding 
generalizes the formula (|22D. Again, the A^-independence of SNq is immediately obvious for 

(3) 

temperatures lower than the restriction temperature , but it actually persists almost 
up to the condensation temperature, as a result of the well-peakedness of the distributions 



V. CONCLUDING REMARKS 

We now come back to our initial question: Can we rely on the grand canonical ensemble 
for Bose-Einstein condensation in isolated traps? As long as we are only interested in the 
usual thermodynamical quantities, the answer is yes. Fig. 5 has shown that the grand 
canonical prediction for the ground state occupation number in a three-dimensional trap is 
almost indistinguishable from the microcanonical one even if the particle number A^ is as 



low as 1000. But in view of the fact that one- dimensional traps are not out of reach |16,38 
it should be kept in mind that these differences are larger for (j = 1 []32 



Another comparison between different ensembles is presented in Fig. 7, which depicts 
the grand canonical (long dashes), canonical (short dashes), and microcanonical (full line) 
specific heat capacity for d = 3 and A^ = 1000. Again, noticeable differences can be found 
only close to the onset of condensation. 

But when it comes to quantum statistical properties, the grand canonical ensemble is 
unsound Adopting the microcanonical spirit and starting from the enumeration of 

microstates, we have computed the microcanonical low-temperature fluctuations SNq for 
harmonically trapped Bose gases. While the case d = 1 could be treated fully analytically, 
leading to eq. (22), the case d = 3 had to rely on the numerical evaluation of the recursion 
relation (p6|). If one could derive an analytical approximation to the numbers Q,^^\n\N) for 
n > N, with a quality similar to that of the Erdos-Lehner approximation ( pT|) to ^l^^^ (n\N) = 
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PN^n), then the case d = 3 could be brought into the same textbook status that the case 
d = 1 aheady has by now. 

It is also worthwhile to point out that very recently Navez et al. have proposed a new 
statistical ensemble, within which the ground state particles act as a reservoir, and exchange 



of particles with the excited-states subsystem without exchange of energy is possible [g9 
The predictions for the fluctuations 5Nq obtained with the help of this so-called "Maxwell's 
Demon ensemble" agree very favourably with our strictly microcanonical results, so that 
this ensemble might constitute a valuable tool for further studies. 

In closing, we reiterate the most important result of the present analysis: while the 
grand canonical ensemble predicts, for ideal Bose gases at low temperatures, ground state 
number fluctuations 5Nq of the order iV, and while usual thermodynamical fluctuations scale 
as \/iV, the low-temperature fluctuations 5Nq for trapped, isolated ideal Bose condensates 
are independent of the total particle number A^. 
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FIGURES 



FIG. 1. Grand canonical ground state fraction Nq/N for a three-dimensional isotropic harmonic 
trap, and gases consisting of 10'^, 10^, and 10^ ideal Bose particles. The larger the particle number, 
the sharper the onset of condensation, and the smaller the shift of the condensation temperature @ 
with respect to Tq = Ig^^ These data result from exact numerical calculations that do not invoke 
the continuum approximation. 

FIG. 2. Grand canonical specific heat capacities for harmonically trapped ideal Bose gases 
{d = 3). The particle numbers are = 10^, 10^, . . . , 10''. The higher the particle number, the 
steeper the drop at the onset of condensation. 

FIG. 3. Microcanonical fluctuations of the ground state occupation number for ideal Bose 
gases in a one-dimensional harmonic trap, for A'^ = 10^, 10^, and 10^ (top to bottom). Full lines: 
data obtained with the help of the saddle point inversion of eq. ([l7|). Dashed: data resulting 
from the Erdos-Lehner formula (]2l]). Tq = Tq^^ = fiujN / {kB^'nN) marks the temperature below 
which the ground state occupation becomes significant |^^. The seeming A^-dependence of the 
low-temperature fluctuations stems from the fact that they have been plotted versus the reduced 
temperature T/T^\ since T^^ depends on N . 

FIG. 4. Microcanonical probability distributions Pe^ {Ncx\n) for finding A^ex out of N = 1000 
isotropically trapped Bose particles excited when the total energy is n • huj. The temperatures 
T/Tq corresponding to the Gaussian-like distributions range from 0.3 to 0.9 (left to right, in 

(3) 

steps of 0.1); the temperature for the rightmost, monotonous distribution is T = 0.95 Tq , which 
is higher than T^^ « 0.93 Tq^\ see eq. (Q). 

FIG. 5. Microcanonical ground state fraction (full line) for A^ = 1000 isotropically trapped 
Bose particles {d = 3) versus reduced temperature T/Tq , compared to the corresponding grand 
canonical data (dashed). 

FIG. 6. Microcanonical fluctuations 6Nq versus temperature, for d = 3 and A^ = 200, 500, and 
1000. The fluctuations are maximal close to the respective condensation points. Below the lowest 
of the condensation points, the fluctuations of the three systems agree perfectly, thus revealing the 
A^- independence of SNq. 

FIG. 7. Comparison of the grand canonical (long dashes), canonical (short dashes), and micro- 
canonical (full line) specific heat capacity for d = 3 and A^ = 1000. 
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